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Abstract 

Diffusion in colloidal suspensions can be very slow due to the cage effect, which 
conhnes each particle within a short radius on one hand, and involves large-scale 
cooperative motions on the other. In search of insight into this cooperativity, 
here the authors develop a formalism to calculate the displacement correlation in 
colloidal systems, mainly in the two-dimensional case. To clarify the idea for it, 
studies are reviewed on cooperativity among the particles in the one-dimensional 
case, i.e. the single-hle diffusion (SFD). As an improvement over the celebrated 
formula by Alexander and Pincus on the mean-square displacement (MSD) in 
SFD, it is shown that the displacement correlation in SFD can be calculated 
from Lagrangian correlation of the particle interval in the one-dimensional case, 
and also that the formula can be extended to higher dimensions. The improved 
formula becomes exact for large systems. By combining the formula with a non¬ 
linear theory for correlation, a correction to the asymptotic law for the MSD in 
SFD is obtained. In the two-dimensional case, the linear theory gives description 
of vortical cooperative motion. 

Special Issue Comments: This article presents methodological insights into math¬ 
ematical treatment of particle cooperativity in colloidal liquids. This article is 
related to the Special Issue articles |m Biophys. Rev. Lett.] about mathemat¬ 
ical approaches to single-file diffusiorP and Fourier-based analysis of quasi-one- 
dimensional systems.^ 
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1 Introduction 


The slowdown of diffusion in systems of densely packed particles, as well as constrained 
dynamics in many other systems described as “glassy”, “jammed”, or “crowded”, has 
attracted attention,!^ not only in connection with the physics of glass transitiorP^ but 
also in the context of anomalous transport in biological cellsP The slowdown is due to 
suppression of free diffusion, which implies that the particles can diffuse only in some 
cooperative way with space-time correlations 

For concreteness, we consider a colloidal system consisting of spherical Brownian 
particles in the Ud-dimensional space, with the representative diameter a. We denote 
the position vector of the i-th particle with ri{t); there are N particles so that i = 
1, 2,..., TV. The particles are governed by the Langevin equation. 


mri 


Vi = -fi-ti - ^ X! 


j<k 


( 1 . 1 ) 


where m is the mass of the particle, Vjk is the interaction potential between the j-th 
and fc-th particles, and is the random force characterized by the variance 

im ® - i')l- (1-2) 

The drag coefficient /i, which may depend on the configuration of the particles in 
general, is regarded here as constant for simplicity. The system is assumed to fill a 
periodic box of the linear dimension L and to be uniform statistically, unless specified 
otherwise, and we consider the thermodynamic limit, L —>■ oo, keeping the mean 
density po = N/L'^’^ fixed. 

To study slowdown of diffusion, let us start with denoting the displacement of the 
j-th particle with Rj(t, s) = rj{t) — rj{s), where s < t; if the system is statistically 
steady, it suffices to consider Rj(t) = rj{t) — rj(0). The mean-square displacement 
(MSB) is given by (R^) = ([Rj(t)]^), with ( • ) denoting the ensemble average in 
regard to the random force, as in Eq. and the subscript j is dropped because 

the system is uniform. In the case of free diffusion, the MSB grows basically as 

(R^) = 2n^Dt (free diffusion), (1.3) 

except for very short inertial timescale, where D = k^T/p is the diffusion constant. 

Suppression of free diffusion in densely packed systems changes the behavior of 
the MSB significantly. Typical behavior of the MSB in glassy liquids is depicted 
in Fig. UKa). After a short period of nearly free diffusion, a plateau appears in the 
MSB, as the particle is trapped in a cage formed by its neighbors. Later the particle 
escapes from the cage and the MSB slowly diverges to the infinity. All this behavior 
reflects different regimes of relaxation in glassy liquids, referred to as the a and /3 
relaxationsThe slowest process is the a relaxation, considered to be the structural 
relaxation that enables the longtime growth of the MSB to infinity. Note that this is 
not a finite-size effect, as we are considering the limit of L —^ -l-oo. The second slowest 
process, corresponding to the period in which the MSB crosses the plateau, is termed 
as the /3 relaxation. 

The slowness of the a and /3 relaxations is expected to have its origin in coopera¬ 
tive motions of the particles, correlated both in space and time. Numerical simulations 
have revealed spatiotemporally nonuniform relaxation referred to as dynamical hetero- 
geneity,l^l^ as well as vortical motion^SHUl and string motions.l^^ In order to capture 
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Figure 1: Behavior of the MSD, showing two-step relaxation. The time t and the MSD 
are nondimensionalized with /D and respectively, (a) A three-dimensional sys¬ 
tem, densely packed with particles subject to Newtonian dynamics with thermostat.^ 
(b) SFD with overtaking allowed [Vmax/^BF = 1, 2, 3,..., up to 8; see the text below 
Eq. (12.1|) ]. The insets show the MSD compensated by t and '/i. 


such cooperative relaxation processes, various space-time correlations have been de¬ 
vised. A popular idea is to define four-point space-time correlations by analogy with 
spin glasseswhile there are different ideas that involve particle tracking, such as 
the bond breaking correlation.1^1^121] Combining the idea of continuum-based particle 
trackin^^ with that of four-point space-time correlations,^ we find that displacement 
correlatiori^^^ can be quite useful. 

The present article aims to give some insights into space-time correlations charac¬ 
terizing cooperativity in glassy liquids, by studying a simpler model of caged dynamics. 
By “a simpler model” we mean the one-dimensional (ID) version of Eq. (II.Ill , whose 
behavior is known as single-file diffusion (SFD). An example of SFD with large but 
finite interaction potential is shown in Fig. [TKb). Here the MSD exhibits two-step be¬ 
havior, similar to the one shown in Fig. [TJa), though there is a difference that the “/3 
relaxation” in standard SFD takes the form of a subdiffusive regime with (i?^) oc \/t 
and not a platearj^- For the limiting case of infinite interaction potential, the a re¬ 
laxation disappears so that the ID cages become eternal. If L is finite, eventually the 
correlated motion dominates the whole system and the MSD grows again in proportion 
to t, though much more slowly than in the free diffusion.l^SHlll 

The article is organized as follows. Basics of SFD are reviewed in Sec. [U with 
emphasis on cooperativity rather than anomalous diffusion. After a survey of certain 
approaches to SDF in Sec. O we introduce in Sec. 0] what we call the Alexander- 
Pincus (AP) formula,l^^l^ to calculate the displacement correlation as well as the 
MSD as its special case. Improvements over the original version of the formulgP^ are 
discussed and the formula is extended to the two-dimensional (2D) case, which is the 
main result of the present work. Applying the AP formula to the nonlinear dynamics 
of ID colloidal liquids, we improve the calculation of the MSD in Sec. [5] Preliminary 
results in 2D colloidal liquids are discussed in Sec. [51 

“it is known that MSD in some extensions of SFD can be slower than see Subsec. 2.1. 
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2 Basic Properties of Single-File Diffusion 


2.1 Subdiffusive behavior 

SFD is an old problenP^II^ but remains an active topic. In particular, there are a 
number of works on SFD motivated by interest in glassy dynamics jU P^Ps JIMllll SFD 
is also related to a study on polymeric entanglemenlPSl and the idea of a figure-eight 
model for kinetic arrest.^ 

Here we consider the Brownian SFD, governed by the Langevin equation, 

mX, = ^ V V(Xk - X,) + (2.1) 

3<k 

which is the ID version of Eq. CH), with Xi = Xi{t) denoting the position of the 
i-th particle. The interparticle potential V (r) defines the particle diameter cr, and 
V^nax = max V (r) can be either finite or infinite. 

Denoting the displacement of the j-th particle with Rj{t) = Xj{t) — ^^^(0), let us 
review the behavior of the MSD, in the overdamped regime {t ^ m/y). Setting 

V{r) = 0 results in free diffusion subject to Eq. (II.dL while the case with I4iax = +oo, 
in which overtaking is forbidden, is subdiffusive. The long-time behavior of the MSD 
for I4iax = +00 is described by the Hahn-Karger-Kollmann (HKK) law,l^l^ 


where S = S{k ^ 0) is the long-wave limiting value of the static structure factor 




(2.3) 




and D’^ = D’^(k —>■ 0) is that of the collective diffusion coefficient, D'^{k), given by 

= D/S in the absence of hydrodynamic interaction.l^l^ol Large but finite values 
of Vmax, as well as finite channel width in quasi-lD systems, allow overt akingj^SSU 
which introduces the “a relaxation” resulting in the two-step behavior in Eig. [TJ De¬ 
tails of the transitional behavior in quasi-ID systems depend on the softness of the 
confinementThe effect of the confinement has been studied through the cross- 
sectional density profil^^ and also through calculation of the free energy barrier from 
the configurational integral.l^l^ 

It may be worth noting that there are several extensions of SED to which the 
HKK law (12.21) does not apply. The HKK law is premised on statistical uniformity 
of the system, absence of long-range interactions, and validity of description by the 
Langevin equation with constant /i and uncorrelated random forcing. Without 
these premises, the HKK law may well be modified. Elomenbom and TalonP^ studied 
the case in which the initial distribution of the particles is inhomogeneous, in such a 
way that one tagged particle is located at Xq = 0 and other particles are arranged 
around it according to the relation X±j oc with 0 < a < 1, so that the 

density decreases as the distance from the tagged particle increases. In this case, 
([i?o(0]^) grows in proportion to which interpolates between the standard 

HKK law (12.211 and Eq. (11.31) for free Brownian particles.l^ The analysis was extended 
to the cases of a heterogeneous single file with distributed diffusion coefhcienfP^ and 
“anomalous” single files whose elementary process involves power-law distribution of 
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Figure 2: Experimental setup of models with granular particles (iron balls with a = 
3 mm) to simulate diffusion of ions through a channel; adopted from Fig. 8 of Ref. [30! 
with modification. The gap length is minimal in (a), while the gap is elongated in (b) 
to form a long pore. 


waiting timep^®^the MSD in these cases can be slower than and even becomes 
as slow as (logt)^ in extreme cases.l^ Some other mechanisms, such as blocking at a 
junctioiP^or spatially correlated noise,can make the diffusion so slow that a plateau 
may appear in the logarithmic plot of the MSD. It is also suspected that the HKK 
law (|2j2|) may break down in regard to longtime behavior of extremely dense single 
files.^l^ In the following discussion of SFD, however, we will mainly focus on the 
standard case in which the HKK law is valid. 

For later convenience, we introduce the length scale £oj expressing the mean dis¬ 
tance between neighboring particles, by po = With p introduced, we 

also define D, D /(.^ and D% / £q, so that 1/D% gives a timescale at which pDH 

becomes comparable to p. For ID systems of rigid spheres, we have £o = 1/po = T/iV 
and D%{k) ~ Po-^[l -f (2po/^) sinafej. 


2.2 Cooperativity in single file 

While many studies on SFD now focus on the anomalous diffusion, it should be empha¬ 
sized that another aspect of SFD was noticed as early as in 1955, when Hodgkin and 
Keyne^^ reported their experiments on potassium ion permeability across the mem¬ 
branes of giant axons from Sepia officinalis. They measured the K+ fluxes through 
the membrane in opposite directions, Qout-s-in and Qin->-out, and compared their ratio 
with Ussing’s equatioiP^ (reminiscent of the fluctuation theorenP^ : 


^out—vin 

Qin^out “ ‘"''P ’ 


(2.4) 


where A$ = q{Eout — Ein) — fcBFlog(pin/pout) is the electrochemical potential dif¬ 
ference, with q denoting the charge of K+; the electric potential and the ion density 
on the side a are denoted with Ea and pa, respectively. Independent diffusion of ions 
should give n = 1, while the experimental value was n = 2.5 > 1, indicating some kind 
of cooperativity. “In a stroke of genius” (if we borrow the words by Hill^^, forty 
years before crystallographic determination of ion channel structure, they said that 
the experimental result can be explained by assuming that K"*" ions tend to move in 
narrow channels. 

To demonstrate this idea, Hodgkin and Keyne^^ performed a mechanical experi¬ 
ment with granular models illustrated in Fig.|3| Causing random motion of particles by 
shaking the system, they observed fluxes between the chambers 1 and 2 and compared 
their ratio with Eq. (El), which reduces to = {pxjppp for electrically 
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Figure 3: Cooperative motion in ID systems, (a) Schematic description of coop¬ 
erative motion due to a diffusing vacancy, (b) Numerically calculated two-particle 
displacement correlation {RiRj) in SFD for different values of f4iax, with N = 10^, 
po = N/L = 0.2(7“^, and the time interval t = 200 a^/D; cited from Ref. HU The 
analytical solution is given by Eq. (EZl). 


neutral particles. The experiment with a short gap in Fig. HJa) gave the ratio close 
to pijP 2 = 100/50 = 2, implying nearly independent diffusion (n ~ 1). Contrastively, 
for a longer gap in Fig. (Hb), Qi _>2 was 18 times greater than Q 2 ->.i, indicating that 
about four particles were cooperating (n = log 2 18 ~ 4.2). 

In regard to ion transport across membranes, researchers have also investigated 
possibilities other than the long pore, assuming cooperativity on the side of the trans¬ 
porting proteins and not among the transported ions, which, however, turned out to 
be unsuccessful.^ The “knock-on” mechanism among the K+ ions is already coopera¬ 
tive, and leads to the high sensitivity of the flux ratio to A<i>. In this case, the system 
size is finite, and the correlation length can span the entire system. 

In situations with L —^ -l-oo, including the case in which the HKK law (El holds, 
the cooperativity exhibits itself through slow diffusion. Cooperativity is usually de¬ 
tected by comparing two configurations at different times, as is illustrated in Fig.^a). 
Here, the leftward diffusion of the vacancy induces rightward cooperative displacement 
of particles. This cooperative motion can be quantified by calculating displacement 
correlationj^l^Sl given by {RiRj) in the ID case. 

A numerical example of displacement correlatiorPI is shown in Fig. Eb). The 
cases with I4iax = k^T and I^nax = are compared. After sufficiently long 

equilibration, the particles were renumbered consecutively at t = 0, and the values 
of {RiRj) were calculated for t = 200a^/D. In the case of relatively low barrier 
(Imax = k^T), the displacements of the neighboring particles are nearly uncorrelated, 
indicating that the particles are diffusing almost freely. Contrastively, in the case 
of higher barrier (Idnax = S/cbT), positive correlation is observed within some finite 
distance. This correlation demonstrates dynamical structure behind the slow diffusion. 
The dynamical correlation length is diffusive and given by \{t) = 2\/DH. 

3 Theories of SFD: extensible to higher dimensions? 

There are various theoretical approaches to SFD. Some techniques, such as the Jepsen 
line mHUl highly depend on the ID geometry, while other techniques may be extensible 
to three-dimensional (3D) colloidal liquids. We look for the latter kind of approaches. 
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requiring in addition that the theory should incorporate cooperativity. 

For the sake of simplicity, in the present section, sometimes we ignore the difference 
between D and — D/S. This is strictly valid only in the limit of point particles 
(ct <C £o = L/N), but it is often possible to take the effect of finite a into account by 
reviving S in an obvious way. 


3.1 Phenomenological treatments of cooperative diffnsion 

Let us start with reviewing the idea of RallisorP^ that relates the slow diffusion to 
cooperativity. If n Brownian particles are strongly interacting and moving together, 
the diffusion coefficient for their center of mass is D/n, in the sense that 


_d 

dt 




2D 

n 


(3.1) 


To apply Eq. (EB to the collective dynamics in SFD, RallisorP^ proposed to replace 
n in the denominator with Af{X) = 1 + poX, which is the number of particles within 
the dynamical correlation length A = A(t), so that the MSD is given by 


^ = ^dirDH. (3.2) 

Upon integration, Eq. (EH) yields 

= —J— -^log (l + poVdirDH^ ; (3.3) 

po V tt 7rpg V / 


it gives free diffusion for small t, and subdiffusion for large t with a logarithmic cor¬ 
rection term. 

It is interesting to test Rallison’s idea in the case of a single file with non-uniform 
initial condition,given by Xo(0) = 0 and X±j(0) oc Since A/’(A) oc A^““ oc 

t(l-a)/2 

in this case, Eq. (13.211 readily yields (i?^) oc xhis result was derived 

by Elomenbom and Taloni^ with an analogous but slightly different approach, in 
which A/’(|.^ was considered instead of A/’(A). 

RallisorPS applied this idea to the 3D colloidal liquid and obtained a critical volume 
fraction, (/c = 0.64 (with ip related to po by (p = |7r(^)^po)j above which the MSD 
cannot grow to infinity. This is greater than the experimental value of glass transition 
point ((^giass = 0.56 according to Pusey and van MegeiP^, and rather near to the 
random close packing density. The discrepancy casts doubt on quantitative validity 
of the theory, but the idea of relating the slow diffusion to the number of particles in 
cooperative motion remains quite suggestive. 

Another instructive idea is found in the lattice SFD by van Beijeren et On 

the basis of the picture of migrating vacancies, i.e. the discrete version of Fig. Ha), 
they calculated velocity autocorrelation of the tagged particle as 

{u{t)u{0)) ~ (3.4) 


for large t, where u denotes the velocity and A„ is a constant that depends on the 
hopping rate, concentration, and lattice constant. With the aid of the relation 


d^ (R^) 
dD 


= 2{u{t)u{0 )), 


(3.5) 
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Eq. (13.41) readily leads to the asymptotic law, oc y/i. 

The velocity autocorrelation in Eq. dSH) has a negative longtime tail: the particles 
are always pushed back by its neighbors. The negative tail can be derived also on the 
basis of continuum description,!^ as will be reviewed later in Eq. (13.1511 . 

Within the framework of lattice dynamics, the idea of migrating vacancies can be 
applied also to geometries other than the ID chain, such as double-file diffusioiP^ and 
2D regular lattices.!^ The idea of clarifying slow dynamics by moving defects has been 
incarnated also in free-volume kinetic models of glasses and granular matters.!^!!^ An 
interesting question is whether this approach can be generalized to systems without 
lattice, such as the one shown in Fig. [3ja). The answer is affirmative in our opinion, 
as will be explained in the latter half of this article. 

The theory of RallisorP^ and that of van Beijeren et were immediately applied 
to SFD on a ring with finite L. The longtime behavior in this case is given by Eq. (EH) 
with n = N. 

3.2 Density correlation in the Fourier space 

As the dynamical correlation length in Brownian SFD is diffusive, it is natural to seek 
indication of collective motions in the density field subject to a stochastic diffusion 
equation. Here we review several ideas for relating the motion of a tagged particle to 
the density fluctuations of the medium, including the celebrated theory by Alexander 
and Pincus.!^ 

Let us define the (microscopic) density field by 

p = p{x,t) = '^Pi, pi = pi{x,t) = 5{x - Xi{t)) (3.6) 

i 

where pi is the single-body density of the f-th particle. Evidently, p satisfies the 
continuity equation 

dtp + d^Q = 0, (3.7) 

with the flux Q = pu given by 

Q = Q{x,t) ='^p^{x,t)Xi{t). (3.8) 

i 

Anticipating a linear equation for p{x, t) as a starting point and assuming that the 
system is translationally invariant, we introduce the Fourier representation: 

P{k,t) = ^ J [p{x,t) - po]dx, p{x,t) = po + '^p{k,t)e~^'^^, (3.9) 

k 

with k/Ak G Z where Afc = 27r/L. 

The two-time correlation of p is referred to as the intermediate scattering function 
and essentially equivalent to the dynamical structure factor,!^ as 

= (^7^0); (3.10) 

the self part of F gives the Fourier representation of the diffusion propagator, 
Fs{k,t)=L^ {p,{k,t)p,{-k,0)) = 

8 


(3.11) 




where pj{k,t) = L~^ exp[ikXj{t)] is the Fourier component of pj{x,t). 

Now let us proceed to the dynamics of p. If the particles are subject to the Langevin 
equation m and attention is focused on the timescale longer than m/p,, then p is 
governed by a nonlinear stochastic equation referred to as the Dean-Kawasaki equa- 
tion.iSSHlll Here we discuss the linear dynamics only, deferring the nonlinear case until 
Sec. [5] The linear stochastic equation for p, sometimes referred to as the diffusion-noise 

equation, EE3 reads 

dtp{x, t) = Ddlp + fp{x, t), (3.12) 

with the statistics of the noise term given by 

{fp{x,t)fp{x',t')) = 2Dda:do,'pix,t)S{x - x')6(t - t') 

~ 2pQDdxdx'5{x — x')5{t — t'). (3.13) 


The linearized Dean-Kawasaki equation (13.121) is readily solved in the Fourier rep¬ 
resentation. With this solution, Taloni and LomhollP^ calculated the velocity auto¬ 
correlation, {uj{t)uj{0)) with Uj = Xj, which is related to the MSD by Eq. (13.511 . The 
single-file condition implies that the velocity of the tagged particle equals that of the 
medium around it (for the timescale longer than l/D,), so that 

u,{t) = u{X,{t),t) c, u(X,(0),t) = (3.14) 

The replacement of Xj{t) by ^j(O) is justified by the relative smallness of the dis¬ 
placement in comparison to the dynamical correlation length. Using the continuity 
equation (13.71) and the space-time Fourier transform of the noise term, they founcP^ 


(u(t)u(O)) 


1 nr 

4po V ’ 


(3.15) 


reproducing the same negative longtime tail as in Eq. (El- We note that the negative 
tail in the velocity autocorrelation is a manifestation of the cage effect, found in 2D and 
3D colloidal systems as well ,1^1^ and differs from the positive longtime tail reported 
by Alder and Wainwright.l^ 

Taloni and LomholtP^ also succeeded in interpolating the longtime collective dy¬ 
namics in SFD and the short-time single-particle behavior, describing the velocity of 
the particle with a generalized Langevin equation 


du 

'dt 


= -[ K{t - t')u{t')dt'+ fKit), 
Jo 


(3.16) 


where /k is a time-correlated noise related to the kernel K by the relation 


{fKit)fK{t')) = mkBTK{\t - t'l). (3.17) 


The longtime behavior of K can be determined by Fourier analysis of the collective 
dynamics in Eq. (13.121) or its chain-dynamical equivalentFurthermore, the 
short-time behavior should reduce to that of a free Brownian particle. By incorporating 
these two limiting cases into K, Taloni and LomholtP^ obtained 


(i?^) 



g4D.t 

-erfc 

2 


(2 



1 

2 


(3.18) 
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Contrastively, the concise theory by Alexander and Pincu^^ avoids directly dealing 
with the velocity u and relates the MSD directly to the density fluctuation, by the 
integral 

(3.19) 

We refer to Eq. (111.191) and its extensions as the Alexander-Pincus (AP) formula. 

The original derivation of Eq. (I11.19|) is traced as follows: Let hj = hj{t) denote 
the displacement of the j-th particle from a configuration in which the particles are 
uniformly distributed, so that hj(t) = Xj(t) — jio- Rewriting p{k,t) with hj as 

p(fc,t) =Y.Pj{k,t) = At)] (fc ^ 0) (3.20) 

j 3 

and introducing h{k,t) = so that hj(t) = , one ex¬ 

pands p{k,t) in Eq. (13.2011 in power series of khj{t), to find 

p{k,t) = [l + ikhj{t) + 0{k‘^h'^)] ~ \kp[)h{k,t)- (3-21) 

j 

By noticing Rj{t) = hj{t) — hj{0) and using Eq. (I3.21|) . the MSD is calculated as 

(i?l) = 2{{[h,iof) - {h,{t)h,m} (3.22) 

k 

then, replacing the summation with an integral and substituting F{k,t) oc *, we 

arrive at Eq. (13.191) . from which (i?^) oc is readily obtained. 

It should be noted that the above derivation includes a step requiring attention. 
Since h is not small, the validity of the expansion in Eq. (I3.21|) poses a delicate issue. 
However, this difficulty is avoided by using a convected coordinate system, as will be 
shown in Sec. |H 

3.3 Mode-coupling theory 

Conceptually, the generalized Langevin equation p.l6l) could be regarded as projection 
of the A'-body equation onto the single-particle motion. However, this projection 
seems to be difficult to perform systematically in the real space (as oppose to the 
Fourier space) It is then interesting to notice the case in which the Mori-Zwanzig 
projection formalism has been applied to 3D dynamics of glassy liquids, known by the 
name of mode-coupling theory (MCT).li^l^l^ Projection of the equation of motion 
onto the density correlation F yields 

{dt + D^k^)F{k,t) = - [ M{k,t-t')dt>F{k,t')dt', (3.23) 

Jo 

which is then closed by approximating the memory kernel M with a quadratic func¬ 
tional of F. In spite of various limitations, the MCT equation has succeeded in repro¬ 
ducing the two-step relaxation of E'(fc, t) and the growth of the relaxation time, as the 
mean density po increases 

Diffusion of a tagged particle is treated by the MCT equation for Fs, in the form 

[dtF Dk^) Fs{k,t) = - f Ms{k,t-t')dt>Fs{k,t')dt', (3.24) 

Jo 
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with the memory kernel Ms expressed as a bilinear functional of F and F^. In de¬ 
scription of SFD with Eq. (I3.24L however, it is difficult to reproduce the \/t behavior. 
As long as the standard memory kernel is used, MCT cannot predict anomalous dif¬ 
fusion.!^ Mathematical difficulty is located in constructing Ms with proper longtime 
behavior, which is actually more difficult than constructing K in Eq. (I3.16L because 
K can be given before the unknown u is determined, while Ms must be a functional 
of the unknown F and Fg. In modified MCT approaches by EedderJ^ and Abel et 
aZ.,!^ Ms is improved by taking the prohibition of positional exchanges directly into 
account; as a result, the correct t^/^-law is recovered. 

Another possible approach consists of retaining Eq. (13.2311 for F alone, while 
Eq. (|3.24|) for Fg is replaced with another equation for tracers, in which the four- 
point correlation is directly taken into account.!^!^ We will see that an improved 
version of the AP formula is suitable for this purpose.!^ 

3.4 Elastic chains 

The longtime dynamics of SED are known to be equivalent to those of fluctuating 
elastic chains. The linear dynamics are then described by replacing V{Xk — Xj) in 
Eq. (12.11) with the harmonic effective potential, so that the equation of motion reads 

mXi + ^iXi = K (Ai_|_i — 2Xi Xi-i) + (3.25) 

with K denoting the effective spring constant.!^!^ Taking the continuous limit and 
considering the overdamped case (m —>■ 0) for simplicity, we rewrite Eq. (13.2511 in terms 
of with ^ here denoting a continuum analogue of i: 

dth{^, t) = D'dlhi^, t) + Ui, t), (3.26) 

where D' = K/fi and {fhi^, , i')) = — — Note that h can be taken 

equal to the position x itself, as h = x{^, t), or the displacement from the mean position, 
as h = x(^, t)—£oC In any case, the displacement is given by t) = h{^, t) — h{^, 0), 
and the MSD is given by the essentially same integral as Eq. (13.1911 . The Langevin 
equation for elastic chains is known as the Rouse model in the context of polymer 
dynamics.!^ The dynamics of a curve h = on the (^,/i) plane can also be 

interpreted as kinetic roughening of a fluctuating surface,with Eq. (13.261) referred 
to as the Edwards-Wilkinson equation.!^!^ 

While this reduction to elastic chain seems obvious in the case of soft-core interac¬ 
tion, the important point is that also the hardcore interaction can be described by a 
smooth effective potential between the particles, as a result of coarse-graining in time. 
Recently, this point has been discussed both in the ID setupl^ and in the context of 
glassy liquids in higher dimension^I^I^, seemingly by independent groups. The effec¬ 
tive potential, calculated from the partition function that gives the free energy, is of 
entropic origin: the interparticle force is essentially the osmotic pressure of the ideal 
solution, i.e. the pressure of the ideal gas. 

The discrete equation ( 3.25p and its higher-dimensional analogues are amenable to 
normal mode analysiswhich is often performed around some metastable config¬ 
uration. Modes with low frequencies are found to be important. In the case of colloidal 
glass with rid > 2, the presence of the low-frequency modes is related to the character 
of the interaction described by the “bonds”, which can directly sustain compression 
but not necessarily resist shear deformation. In the continuum description, the nor¬ 
mal modes are plane waves. In correspondence to the low-frequency modes in discrete 
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systems, the continuum description of colloidal liquids must respect the distinction 
between the longitudinal and transverse modes. 

The situation can be simpler in systems that are elastic by nature, such as an 
assembly of densely packed softcore particles or a system of colloidal particles trapped 
in a polymeric network. In such models, the presence of finite rigidity could be taken 
for granted, so that simplified description based on the scalar integral 

f \ 

(i?^) oc / - — -d"‘^k (without tensorial treatment) (3.27) 


might be justifiable. Once the description by the linear elastic model is justified, 
it allows analytical calculation of displacement correlation and even some four-point 
correlations such as dynamical susceptibility.^^ Details will be discussed elsewhere. 

The effect of dimensionality on the MSD in elastic models may deserve some com¬ 
ments. The ID AP formula is free from ultraviolet divergence, while t must be finite 
to avoid infrared divergence of the MSD. In the 2D and 3D cases, a ultraviolet cutoff 
must be introduced. Denoting the cutoff wavenumber with we calculate 

the MSD in the 2D case from Eq. (13.271) as 


(R^) oc log (^1 + klDli^ , 


(3.28) 


and for na = 3 we find 


(R2) oc 1 - 



(3.29) 


The large t behavior of Eq. (13.281) differs from that of Eq. (13.291) : the MSD for na = 3 
converges for t —>■ -|-cx), while that for na = 2 diverges logarithmically, until the effect 
of the finite system size comes into play at t ^ L^/D^. Thus the 2D elastic model 
suffers both infrared and ultraviolet divergences. It can be even argued that the 2D 
harmonic solid is not a solid in the usual sense,!^ at least if one accepts the argument 
that single-file systems are denied to be solid because any particle in the system can 
get arbitrarily far away from its initial position.l^ 


4 Alexander—Pincus Formula in Convected Coordi¬ 

nate System 

We have reviewed several theories of SFD that may be extensible to higher dimensions, 
including the AP formula p.l9l) . which relates the MSD to F{k,t) oc Since 

the original derivation of the formula is of approximate nature, it seems difficult to 
improve the result by considering a nonlinear correction to F. 

Here we redevelop the AP formula, making it exact by changing the variables ap¬ 
propriately.!^^ The key to this approach is the adoption of Lagrangian correlations, 
motivated by theories of the Navier-Stokes turbulence.l^^lf^ While the formula itself 
remains linear with regard to the correlation in the new variables, it can be combined 
with a nonlinear equation of the correlation, so that the result can be improved as we 
will see later in Subsec. The formula is also extended so as to give the displacement 
correlation in rid-dimensional systems. 
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4.1 Label variable and the ID AP formula 


Let us start with reviewing a classic idea of the Lagrangian description in fluid me¬ 
chanics. The idea is to introduce a continuous and invertible mapping 


(Id) 


^ r = 


y(ld) 


(4.1) 


such that its inverse, ^ satisfies the convective equation 


(5t + u-V)|(r,t) =0. 


(4.2) 


This mapping makes it possible to specify a world line as ^(r, t) = const., and therefore 
we refer to ^ as the label variable. In other words, ^ is a continuum analogue of the 
particle numbering; as a matter of fact, the ID label variable is a mere continuous 
interpolation of the consecutive numbering of the particles, but here we introduce the 
label variable in a more abstract way so as to enable extension to higher dimensions 
without postulating regular numbering at all. The moving curvilinear coordinate sys¬ 
tem given by Eq. dun is sometimes referred to as the convected coordinate system^^ 
The description in which the field variables are regarded as functions of (^, t) is termed 
as the Lagrangian description, as opposed to the Eulerian description in which the in¬ 
dependent variables are r and t. 

There are infinitely many choices of label variable in general. In the ID case,l^l^ 
we choose ^ so as to satisfy 


dx^{x,t) = p, dt^{x,t) = -Q, (4.3) 

where p = p{x, t) and Q = Q{x, t) are the density and the flux satisfying the continuity 
equation (HTfl) . It is then readily verified that ^ = ^{x,t) satisfies the ID convective 
equation, which qualifies ^ as a label variable. 

For this particular choice of the label variable, the derivatives of a: = x{^,t) are 

d^x{^, t) = ^ ^ ^ 

The presence of 1/p in Eq. (14.411 motivates us to introduce the “vacancy field” 

so that the relation dtd^x = d^dtx is cast into the form of the vacancy conservation, 

t) = d^u{i, t). (4.6) 

Then, by using '0(|,t) instead of p{x,t), we can derive a more rigorous version of 
the AP formula,!^ in the sense that the approximation in Eq. (13.211) is not required. 
The formula for the MSD reads 

^ C{k,Q)-C{k,t) 

\ / 7riV2 ]_^ k^ 

in the statistically steady case, where 

C{k, t) ^ (-!/>(fc, t)ip{-k, 0)) 


(4.7) 

(4.8) 
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is the Lagrangian correlation of vacancy in Fourier representation: 




The hacek (inverted hat) is used to make a clear distinction between the two Fourier 
transforms exemplified by Eqs. (13.91) and (14.91) . The formula for the MSD in Eq. (j4.7l) 
is a special case of the formula for the displacement correlation, 




ttN"^ 





(4.10) 


where R{^,t) = For ^ Eq. (j4.10ll reduces to Eq. (j4.7|) . 

To begin the derivation of the AP formula, we notice that 


d^R{^, t, s) = £o [fpii, t) - tp{^, s)] 


(4.11) 


where i?(^,t, s) = x{^,t) — x{^,s)-, this is readily verified with Eqs. (14.41) and (14.Sp . 
Integration of Eq. (14.111) in Fourier representation yields 


Rit t, s) = 4 ^ [ip{k, t) - '^(k, s)] 

k 


—iA: 


+ Rcitj s), 


(4.12) 


with Rq denoting the displacement of the center of mass, negligible for N ^ oo. 
Multiplying Eq. (14.121) by its duplicate with (4 k) replaced by (4, —k'), we find 


(i?(4t,s)i?(4,t,s)) 

__ __ _ _ _ _ p-ifcC+ifc'4' 

= ^0 XI XI ( - ^i-k', s)])-- 

k k' 

= -^0 X (“ '^i-k, s)])-p-; (4.13) 

k 


the terms with k ^ k' is shown to vanish if ipik^t) is subject to linear Langevin 
equation whose forcing term is uncorrelated for different wavenumbers [see Eq. (1^ . 
for example], and the same is true for nonlinear cases, in the limit of A^ —^ ck), as 
long as the vertex is sparse in Kraichnan’s sense.^^^ Taking continuum limit as 
Efc(- ■■)^ £/(■■■ )dfc, we obtain 


{R{^,t,s)R{C,t,s)) 

Ri r+°° \C{k,t,t) + C{k,s,s) 
~^J-oo [ 2 


C{k, t, s) 


P 


dfc 


(4.14) 


where C{k,t,s) {N/L'^) (^'tp{k,t)'ij){—k, s)). In the statistically steady case, we can 
choose s = 0 without loss of generality and, setting C{k,t,0) = C{k,t), which also 
implies that C{k,t,t) = C{k,s,s) = C'(fc,0), we arrive at Eq. (I4.10F 

The AP formula (14.101) has the advantage of being able to provide the displace¬ 
ment correlation, which is a four-point space-time correlation, in terms of two-body 
Lagrangian correlation. This is contrastive to the approach via Fg in Eq. (13.111) , which 
requires calculation of a four-body correlation. 
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4.2 AP formula in higher dimensions 

The AP formula may seem to be easily extensible to higher-dimensional cases, simply 
with the rid-dimensional wavenumber integral as in Eq. (jd.27ll . This form, however, 
should not be taken literally but only as a simplified treatment in which the deforma¬ 
tion is restricted to one of the modes of elastic wavesP^ As was noted in Subsec. [311 
the description of colloidal liquids requires distinction between the longitudinal and 
transverse modes. Here we demonstrate how to extend the AP formula respecting this 
distinction, which results in a tensorial formula. 

For concreteness, let us focus on the 2D case (rid = 2). Our target is the displace¬ 
ment correlation tensor, 


(R(|,<,s) 0 R(|',t,s)) 


(RiRi) {R1R2) 

{R2R1) {R2R2) 


(4.15) 


where i?i and R 2 are the Cartesian components of the displacement R, with the 
arguments omitted when obvious. The label variable ^ = (^, 77 ) is defined through 
the 2D extension of Eq. (HD) . later shown as Eq. (14.201) . so that = ^{Yi{t),t) is 
independent of t. Conceptually, (■oE.i is supposed to express the mean position of the 
i-th particle (for some appropriate timescale shorter than the a relaxation tim^^. 

The ambiguity due to non-uniqueness of ^ can be removed by introducing 


X{d,t,s) 



(5^ (re - ta - d) Ra 0 Rb d^rAd^re 


(4.16) 


where ta = r^A, s), Ra = ^(^a, t) — r(^A, s), and tb and Rb are defined analogously. 
Equation (I4.16P gives the two-particle displacement correlation as a function of the 
initial relative position vector d. Approximately we have 


X(d, t, s) ~ (R(|, t, s) 0 R(|', t, s)) (4.17) 

with d ~ .^o(^ ~ € 0 ) and a correction to this approximate expression can be obtained 
in terms of triple correlations (see Subsec. V-A. in Ref. 1^. 

Recalling that our construction of the ID label variable through Eq. (1131) is closely 
connected to the ID continuity equation (EH), we relate the rid-dimensional label 
variable to the rid-dimensional continuity equation. 


dtp -l- V • Q — 0, 


(4.18) 


in a similar way. The (particle-scale) density p = plv, t) = Pi and the corresponding 
flux Q = Q(r, t) = pu are given by straightforward generalization of Eqs. (13.61) and 
(13.81) : for rid = 2, we have pi = (5^(r - Yi{t)) and Pi(r, t)ri(t). 

To generalize Eq. (14.31) to rid-dimensional cases, we write it in the form 


{p, Q) 


Go 

ei da;^ 


{dx^, -dtO- 


(4.19) 


The above form is suggestive and motivates us to assume that ^ is related to (p, Q) 
by the bilinear equation 



60 dt^ dtp 

(p,Q) = 

dx^ dxP 


02 dyS, dyP 


C ^x V 


dy^ dyP 

dyT] 

5 

dt^ dtp 


dt^ dtp \ 
dx^ dxP J 


(4.20) 
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for rid = 2, and a trilinear equation in an analogous form for the 3D caseP^ It is then 
verified that ^ satisfies the convective equation (14.21) . if it solves Eq. (14.201) . Within 
linear approximation in regard to p — po and Q, a solution 

(e, v) = - 4 f Q(r, t) di + 4V0, = p(r, o) - po (4.21) 

«0 Jo 

can be constructed, where o is an arbitrary “initial” time. 

In the spirit of Eqs. 631) and 631) where t/j = ' 0 (^ 1 ^) is introduced on the basis of 
dx/d^ = 1/p, we consider the deformation gradient tensor and introduce ( 4 / 1 ,'I' 2 ) on 
the basis of its diagonal components, as 


d^X drjX 
diV drfV 



* 


* 


(4.22) 


where the off-diagonal components are omitted and replaced with an asterisk. 

It follows from Eq. (14.221) that s) = £0 ['I'i(£,t) — 'I'i(^,s)], with s < t, 

and an analogous expression for driR 2 is also readily available. Introducing the Fourier 
representation as 


■^i(et) 

,«' 2 (l,t). 


E 


5'i(k,t) 

4'2(k,<) 


„-ik€ 


(4.23) 


with k = (/ci, ^ 2)5 which is a straightforward generalization of Eq. (14.9L we find 


k 


g-ik€ 

—iki 


+ Xcit, s); 


(4.24) 


subsequently, with the fluctuation of the center of mass Rq = (^gWg) neglected for 
N ^ 00 , Eq. (14.241) yields 

(Ri(|,t,s)Ri(|',t,s)) 

= ^oEE(['I^i(k,t)-4/i(k,s)] [4/i(-k',t)-vI/i(-k',s)])^--. (4.25) 

k k' 

Provided that the same condition as was mentioned immediately after Eq. (14.131) is 
satisfied, the terms with k ^ k' are seen to vanish as in the ID case. Then, taking 
the continuum limit, we obtain an expression that gives (i?ii?i) in terms of correlation 
of 'ki. Other components of the displacement correlation tensor are calculated in a 
similar way. Defining 

Co,p{k,t,s) = ^{^o.{Kt)^p{-Ks)) (a,^G {1,2}) (4.26) 

and assuming Cap to be real and symmetric so that C' 2 i(±k, t, s) = C'i 2 (k, t, s), we 
find 


{Ra{$,t,s)Rp{$' ,t,s)) 


2Tr^N 


Capikj Sj s) Cap{k-, t, t) 


- Cap{k,t,s) 


3 -ik.(C-«') 


-dfcid/c2, (4.27) 


where we can choose ^' = 0 without loss of generality. 
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For later convenience, let us rewrite the 2D AP formula in terms of the dilatational 
and rotational modes of deformation. We define 


'(/'d(k,t) = ^'i(k,t) + ^2(k,t), (4.28a) 

^/’r(k,t) = ■^'f'i(k,t) - ■^^'2(k,t), (4.28b) 

/C2 k,i 

so as to have the Fourier representation of the dilatational and rotational modes: 


d^x + d^y = 4 


2 + X!^d(k,t)f 


-ik€ 


= 4 ^V’r(k,t)e (4.29) 

k 


The correlations of these modes are denoted by 

N 

Cab(k, = JJ ('*/'a(k, t)'!/'b(-k, s)) (4.30) 

with a, 6 G {d, r} and s <t. Obviously, Cap in Eq. (14.261) can be expressed as a linear 
combination of Cat in Eq. (14.3011 . and vice versa. For simplicity, let us assume that 
Crd and Cdr vanish identically for some reason (such as the parity), and denote Cdd 
by Cd and C^r by Cr- Then the relation between Cap and Ca reads 


Cii 412 
421 C22 


1 

io 


■ kt 

A/2 


U‘2u2 

^1^2 


4d 




1 -1 

-1 1 


4r; 


(4.31) 


note that Cap is symmetric. Substituting Cap in Eq. (14.311) into Eqs. (14.2711 and 
choosing = 0 without loss of generality, we find 


(i?i {i, t, s)Ri (0, t, s)) {Ri (^, t, s)i?2(0, t, s)) 
(i?2(l, t, s)i?i(0, t, s)) (i?2(l, t, s)R 2 { 0 , t, s)) 


L» // 

■4d(k, s, s) -h 4d(k,t,t) 

27r2A JJ 

2 

L> ff 

■4r(k, s, s) -h 4r(k,t,t) 

27t^N JJ 

2 


- 4d(k,t, s) 
4r(k,t, s) 


■ kf 

—4^1 


fci4 

—44 

kl 


g-ik-S 
,-ikC 


k4 


d/cid4 
dA:id4- 


(4.32) 


Equation (14.3211 is our main result in this article. 

As we will see in Sec. [51 Eq. (I4.32|l is simplified in the isotropic case, in which 
4d and 4r are independent of the directions of k so that Ca = Ca{k,t,s). In this 
case, with the approximation d ~ 4(^ — €0) Ike displacement correlation tensor is 
expressible in terms of two functions Ay = A||(4t) and Aj_ = X±{^,t): 


X = A||(J/4,t)^ + Aj.(d744) ( 1 -^) • 


(4.33) 


The functions Ay and Aj_ denote the longitudinal and transverse displacement corre¬ 
lations, respectively. 


5 Displacement Correlation in ID Colloidal Liquid 

To demonstrate how to apply the AP formula to colloidal systems, here we begin by 
illustrating the ID problem. Starting from the Dean-Kawasaki equation, we calculate 
the Lagrangian correlation C{k,t) and substitute it into the ID AP formula. This 
procedure will be a prototype for the 2D calculation to be discussed in Sec. [SI 
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5.1 Formulation 


The dynamics of Brownian particles in the overdamped regime can be expressed as a 
stochastic equation for the density field, known as the Dean-Kawasaki equationl^^Hlll 
Utilizing the continuity equation (lT7l) that relates the density field p{x, t) to the flux 
Q{x,t), we write the ID Dean-Kawasaki equation in the form 


Q = -D 



j 


(5.1) 


where U = U[p]{x) = f Vf,a{x — x')p{x')dx' describes the interaction of the particle, 
with VeS denoting the effective potential based on the direct correlation function. 

Subsequently, we rewrite Eq. (EU with the label variable ^ introduced through 
Eq. (14.31) . The field variable is ijj = and the corresponding flux, according 

to Eq. (14.61) . is u = Dividing both sides of Eq. (15.11) by p and rewriting the 

differentials with the chain rule, dx = pd^, we obtain 


u = —D 


= -D 




(5.2) 


where /l(C, t) = 9^ E* ^ “ “*) /*(^) and 

Switching over to the Eourier representation given in Eq. (HU), we rewrite the 
combination of Eqs. (ira and (1^ as 


dt'4}(k,t) =-Dlk^'^(k,t)+ ^ Vl'^'ipi-p,t)ip{-q,t) + pofi,{k,t), (5.3) 

p+q+k=0 

with terms of 0{'tp^) dropped. The concrete form of the vertex is given in Refs. [55] 
and 1231 The statistics for the thermal forcing are given by 

9D 

pl{fL{k,t)fL{-k',t')) = -j^k‘^5kk'5{t-t'). (5.4) 


5.2 One-dimensional linear theory 

As an input into the AP formula, let us calculate C{k,t) in Eq. (14.81) on the basis of 
linear approximation to Eq. dOj). We assume that the system is in equilibrium, so 
that the initial condition is not important. Then the linearized equation yields 

ip{k,t) = f dt'e~^°'"^^*~*'^pofhik,t’), (5.5) 

J —OO 

from which we obtain, with the aid of Eq. (El, 

C{k, t) = S' = S{k) ~ S(0). (5.6) 

Using the the AP formula (14.101) . we convert the vacancy correlation C into the 
displacement correlation, ,t)). The result is shown to be expressible in 

terms of the similarity variable 9 = £o{^ — with \{t) = 2^/lJH^ as 

{R{^,t)R{C,t)) = (e"®" - V^|6»|erfc|0|) . (5.7) 
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Setting 0 = 0 in Eq. (15.71) reproduces the established HKK law (12.21) . The prediction 
of Eq. (15.71) is compared with numerically computed displacement correlation^^ in 
Fig. El^b), with j — i interpreted as ^ The numerical correlation for the higher 

barrier case (14iax = ^k^T) is consistent with Eq. (15.71) . 

Incidentally, the analysis can be extended to a case of ageing. Suppose that the 
system has been equilibrated, in the presence of an extra repulsive interaction between 
the particles, in the state characterized by the static structure factor S'init(fc)- Sub¬ 
sequently, the extra interaction is switched off at t = 0, so that the system relaxes 
toward a new equilibrium state. Using the linear solution 

ip{k,t) (5.8) 

^0 

and Eq. (El), we find 

C{k, t, s) = - ^ e-Dlk^it+s) 

with 0 < s < t, S'init — 5'init(0) and S ~ S'(O). Note that, as expected, Eq. (15.91) 
reduces to Eq. (15.61) for S'init -t S. 

It is then straightforward to calculate the MSB in this ageing case. By substituting 
Eq. (15.9p into Eq. (|4.I4p with $ = we obtain 


Po 


Sinit - S (jD<={t + s) 
Po \ 




(5.10) 


In studying ageing effects,^ it is more convenient to rewrite the time arguments in 
Eq. (15.101) as s = Tw and t = %, + t' , denoting the “waiting time” with %,■ It 
is then readily shown that {[R{%; + /yfV is a function of t' jT^ only. This 

function is asymptotically independent of Sinit for t' <C Tw, so that the HKK law 
(12.21) is recovered; on the other hand, as t' elapses and exceeds 7^, the effect of Sinit 
reappears. Eor simplicity, let us focus on the case of the initial condition with equally 
spaced particles,!^ which corresponds to Sinit = 0. For this case, Eq. (IS.IOp yields 


([s(rw + t',rw)]^) ^ f 1 {f « Tw) ^ 

{2Slpa)^/D^ (t'>rw), 

showing that, if t' is the same, the ratio of the MSB for long Tw to that for short Tw 
is : I. Thus the present analysis, valid for cr > 0, gives the same factor as was 
predicted by Leibovich and BarkaP^for point particles (cr = 0). 


5.3 Nonlinear analysis of ID colloidal liquid 

Now we return to the statistically steady case and discuss the effect of the nonlinearity 
in Eq. (15.31) which was ignored in Subsec. 15.21 There are two sources of nonlinearity: 
the one originating from the direct contact between particles is unimportant for low 
density so it may suffice to incorporate it via = DjS, but there is another kind of 
nonlinearity that comes from the configurational entropy, whose effect is not negligible 
even in the case of low density. 
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The nonlinearity of the configurational entropy manifests itself as a memory effect 
due to vacancy-vacancy interaction!^ For SFD on a lattice, this memory effect was 
pointed out by van Beijeren et alW^, assuming some phenomenological rules about 
the dynamics of a vacancy cluster, they calculated (i?^) for both short and long t. 
The crossover between the two limiting cases of D^t <C 1 and 1 in the case of 

continuous systems was taken into account, again phenomenologically, in Eq. (1531) by 
RallisorP^ and in Eq. (13.181) by Taloni and Lomholt.l^ 

It is then interesting to ask whether a more systematic treatment of the nonlin¬ 
earity can improve the phenomenological results in Eqs. (13.311 and (13.1811 . Multiplying 
Eq. (15.3|) with ip{—k,t) and developing mode-coupling theory to evaluate the triple 
correlation, we find the Lagrangian correlation C to be governed by 

{dt + Dlk^)C{k,t) = - f Mi^{k,t-t')dt'C{k,t')dt' (5.12) 

Jo 

with the memory function Ml which is a quadratic functional of C. By calculating C 
from Eq. (15.121) and substituting it into the AP formula (14.711 . our groufP^ found 


(i?^> 


2 jlDt y/2 _2 

po V tt 37r 


(5.13) 


In Fig. m the MSD in Eq. (I5.13|l is compared with direct numerical simulation 
(DNS) of interacting Brownian particles. The phenomenological equations (13.31) and 
(13.181) are also included. Since S = 0.624 in this case differs considerably from unity, 
we have revived S in the equations in an obvious way to make the asymptotic behavior 
consistent with the HKK law (12.211 . 

With ^/Ui/a taken as the horizontal axis, the HKK law (Ell) is expressed as a 
straight line passing through the origin in Fig Ea). Comparing this straight line 
((i?^)jjp,j^) with the result of ((-^^)dns’ bullets), we notice that a negative 

correction is needed. Judging from the graph of (^^)l)ns ~ ('^^)hkk FigHb), the 
correction should become a negative constant for t —>■ oo. 

All the three equations, namely Eqs. (13.31) . (13.181) . and (15.1311 . propose negative 
correction to Eq. (12.21) . Figure H] allows comparing the predictions of these three 
equations with the result of DNS (except for the close vicinity of t = 0 governed by 
the inertial effect). The logarithmic term in Eq. (13.31) is found to give an overcorrection. 
The behavior of Eq. (13.181) seems close to that of DNS for \/lji < Iq = 4cr, but for 
larger t it also becomes overcorrective. Among the three equations, the best estimation 
of the asymptotic behavior is given by Eq. (15.131) . which emphasizes the usefulness of 
the AP formula combined with Lagrangian MCT. 


6 2D Colloidal Liquid 

Finally we have arrived at the stage where analysis of 2D colloidal liquids is within 
our reach. The greatest difference in comparison to the ID case is the presence of the 
transverse mode. 

6.1 Linear analysis of 2D colloidal liqnid 

Here we target on calculation of both the longitudinal and transverse displacement 
correlations, namely Xy and X±_ in Eq. (I4.33L based on linear analysis of the 2D 
Dean-Kawasaki equation. 
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Figure 4: Comparison of the ID MSD calculated by the HKK law (12.21) . Rallison’s 
equation (13.31) , the Taloni-Lomholt solution (13.181) , and the present result in Eq. (I5.13P , 
with the result of where N = 2^® and po = N/L = 0.2bla. (a) The values 

of MSD, nondimensionalized with are plotted against y/Dt/a. The straight line 
shows the HKK law, (R^)hkK’ ™ ^2.21) . (b) The difference between the calculated 

MSD and (i?^)HKK case are plotted in the same way as in (a). 


The approach is basically in parallel with the ID case. We start from the Dean- 
Kawasaki equation, written in the form 


Cl = -D 




3 


( 6 . 1 ) 


where U = U[p]{r) = f 14ff(?’*)p(f)d^f with r* = |r —f|. The thermal noise ij{t) 
is subject to Eq. dni). Changing the independent variables from (r, t) to (4,t) and 
taking the field variable ('I'l, 'I'2) whose conservation law is given as 









we rewrite Eq. (EH) in the form corresponding to Eq. (lOD : 


( 6 . 2 ) 


u = -D 






(6.3) 


where V = (V^)^^ + {\7r])drj is expressible in terms of and 'I'2. 

Linear analysis of Eqs. EH and EH is then performed. In terms of the Fourier 
modes introduced in Eq. (I4.23L the linearized equation reads 


^i(k,t) 

4'2(k;,t) 



kf 

kl 


^i(k,t) 

^2(k,t) 


0-1 


'A(k,t) 

/2(k,t). 


(6.4) 


with the forcing statistics given by 



■(/i(k, t)/i(-k', t')) (/i(k, t)f2(-k', t'))' 

2D, 

'kl O' 

^0 

.(/2(k, t)fi(-k', t')) (f2(k, t)f2(-k', t'))_ 

N 

. 0 kl 


t'), (6.5) 


where the delta-correlation in regard to the wavenumber originates from the random¬ 
ness of the particle configuration; possible correction at shorter length scales is ignored 
for the present. 
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The linear part of Eq. (16.41) has two eigenvalues: fcj + and 0, correspond to 
the dilatational and rotational modes, respectively. The correlations Ca/ 3 ()^,t, s) in 
Eq. (I4.26|) are then calculated, according to Eq. (14.311) . as linear combinations of 
Cd = C'd(k, t, s) and Cr = Cr(k, t, s). These correlations are given as 

Cd(k,t,s) = ^exp s)] , (6.6a) 

on u2 

C,(k,t,s) = -^(s-o), (6.6b) 

where o is a constant of integration, interpretable as the time at which Cr is reset. 
Note that both Cd and Cr are independent of the direction of k. 

The last step is to convert Cd and Cr into the displacement correlation. Using the 
2D AP formula in the form of Eq. (I4.32|) and introducing 9^ = + rf‘)/{AD%t) ~ 

d^/[A(t)]^ with X{t) = 2'/DH, we obtain 


{Ri{$,t,s)Ri{0,t,s)) = 
(i?2(|,t,s)i?2(0,t,s)) = 
(i?i(|,t,s)i?2(0,t,s)) = 


S' 


47rpo 

S 

47rpo 


El (02) 
El (02) 


^2 — e 

^2+ry2 02 

— ^2 -I- jy2 g-e" 
^2 _|_ ^2 02 


S iv e 




Sttpo ?2 +772 02 

which can be rearranged, in terms of </> such that d = {dcoscj), dsincj)), as 


(R(^,t, s) (g) R(0,t, s)) = 


S 


47rpo 


El (02). 


■,-e^ 


02 


02 


■ COS 2(j) 


02 


■ sin 2(j) 


■,-e^ 


(6.7a) 

(6.7b) 

(6.7c) 


■ sin2(/) 


El (02)- —cos2(j) 


( 6 . 8 ) 

By comparing Eq. (16.81) with Eq. (14.331) . on the ground that the expression on the left- 
hand side of Eq. (16.81) equals X according to Eq. (14.171) . we reach the goal of forming 
expressions for and A'j_. For large 0, using the asymptotic form of the exponential 
integral Ei(02), we find 


^11 (eu) 


s 

dTT/Oo 


El (02) 


S 

47rpo 


El (02) 



02 


47rpo 


( 20 - 2 - 0-4 + ...) 




47rpo 


(-0-^ + ---)- 


(6.9a) 


(6.9b) 


6.2 Discussion 

According to Eq. (l6.9aL the longitudinal displacement correlation is positive, while 
the transverse displacement correlation X± in Eq. (I6.9bl) is negative. This is qualita¬ 
tively consistent with the vortical cooperative motion observed by Doliwa and Heuer.li^ 
Numerical verification of Eqs. (EH) will be reported elsewhere.!^ 

By analogy with the ID displacement correlation in Eq. EH whose limiting value 
for 0 —>■ 0 correctly gives the HKK law (j2.2l) of MSD, one may naturally expect that the 
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2D displacement correlation in Eqs. (16.91) will provide the MSD in the limit of 6* —>■ 0. 
Unfortunately, it turns out that the 2D linear theory does not give a quantitatively 
correct value of the MSD, because it is not accurate enough in the short-wave region. 
We may still attempt, however, a qualitative estimation of the MSD in terms of Xy 
in Eq. (I6.9all . Introducing the nondimensionalized cutoff length ,fcut ^ and using 
asymptotic form of the exponential integral for small 0, we may estimate the MSD as 

(R2)^X||(^eut,t)-^0(l)x -^ + ^log^. (6.10) 

Thus we obtain normal diffusion with a logarithmic correction term. This seems to be 
qualitatively consistent with the result of van Beijeren and Kutnei^^ for 2D lattices. 
The logarithmic correction seems to have the same physical origin as the MSD for 2D 
elastic bodies in Eq. (13.281) . 

Nevertheless, we must emphasize that the validity of Eq. (16.101) is rather limited 
and, in particular, it is definitely incorrect for dense colloidal systems such as the one 
depicted in Fig. [3Da). This is not the fault of the AP formula but the limitation of 
the linear theory, as the time needed for the a relaxation is largely underestimated 
by ignoring the nonlinear effects of the interacting particles with finite tr. Evidently, 
a nonlinear theory for Cd and Cr should be developed and used instead of Eqs. (EH). 
Combination of such nonlinear theory and the AP formula will provide a better esti¬ 
mation of the MSD in 2D systems. 

Some insights into the a relaxation may be given by comparative discussion on 
SFD in quasi-lD systems and 2D dense colloidal liquid. The a relaxation in colloidal 
liquid is believed to involve cooperation of numerous particles that occurs in some 
localized yet extended mode, often termed as a cooperatively rearranging region.l^ As 
a drastic simplification of a cooperatively rearranging region in the spirit of mean-field 
approximation, we may replace most of the particles with confining walls and consider 
only two particles explicitly, so that the a relaxation is modeled by the positional 
exchange of two particles in a quasi-lD channel. Let us denote the position vector of 
the particles with ri = (Ai,li) and r2 = (A2,l2), and the potential with 

V = C/(ri, r 2 ) = U(ri 2 ) + C/ex(Ui) + (6.11) 


where U(ri2) is the interaction potential as a function of ri2 = |r2 — ri| and Uex{Yi) is 
the confinement potential for the *-th particle. The barrier potential for the positional 
exchange of Xi and X 2 is given in terms of configurational integral, as 


UiD = UiD(Ar2 — Ai) 


-feeTlog (^JJ e-P^dYidY 2 


const.. 


( 6 . 12 ) 


where /3 = l/^ksT). The “hopping rate” for the positional exchange,!^ denoted with 
l/^hop, is proportional to Let us discuss how Thop depends on Uex, com¬ 

paring two cases: 

Uex(U,) = UsoftiY^) = (6.13a) 


and 


t/ex(U,) = t/hard(U*) 


0 (|r,| < Ly) 

00 {\Yi\>Ly). 


(6.13b) 
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By assuming V (r) to be rigid sphere potential in both cases for the sake of simplicity, 
the integral in Eq. (16.121) can be concretely evaluated. In the case of the soft wall 
in Eq. (l6.13aL we find logrhop oc Ka^/(kBT) in the limit of low temperature; this 
expression of Thop has no singularity as a function of k and a. Contrastively, Thop 
in the hard wall case with Eq. (I6.13bl) has a singularity as a function of Ly and a, 
as it diverges for cr/Ly — >■ 1/2. Thus we find that the cr-dependence of Thop, or the 
dependence of Thop on the strength of the wall confinement, can be sharp or blunt 
depending on whether the wall is rigid or not. This behavior seems to be consistent 
with the numerical findings of Lucena et In the case of 2D and 3D colloidal liquids, 
the corresponding problem is how the timescale of the a relaxation depends on the 
volume fraction </. It is known that this timescale can exhibit an extremely sharp 
dependence on expressible as with B and 0o being some constantsj^l^ 

in comparison, the singularity of Thop for UexiYi) = 17hard(hi) in the quasi-lD channel 
is much weaker, as it diverges only algebraically. This difference suggests subtlety of 
the cooperative dynamics in the a relaxation in colloidal liquids that is difficult to 
capture by the “mean-field” approach discussed here. 

7 Concluding remarks 

In search of insight into cooperativity associated with the cage effect in glassy liq¬ 
uids, we have developed a formalism to calculate the displacement correlation tensor, 
(R(^, t) 0 R(^', t)), in TT-d-dimensional systems. The calculation relies on the label 
variable method, i.e. adoption of the Lagrangian description.l^l^ 

In the ID case, Eq. (gH) gives the MSD as a special case of the displacement 
correlation in Eq. (|4.10|) . These equations improve on the original formula (13.191) in 
that they are asymptotically exact for iV —>■ c» by virtue of the Lagrangian description. 
By combining Eq. 621) with the Lagrangian MCT equation (I5.12|l , we have calculated a 
finite-time correction to the asymptotic HKK law (12.21) . The formula was also applied 
to an ageing SFD and reproduced the effect of the initial condition on the MSD, 
previously obtained by Leibovich and BarkaP^ with a different approach. It will be 
interesting to study various extensions of SFD with the AP formula (14.71) . so that the 
effect of spatially correlated nois^^ may be clarified, for example. 

Our main result is the 2D formula (I4.32|) . which is tensorial. For isotropic systems, 
the displacement correlation tensor reduces to the longitudinal correlation Xy and the 
transverse correlation X±. Linear analysis of the two-dimensional colloidal system 
predicts X± > 0 and Xy < 0, which seems to be qualitatively consistent with the 
vortical cooperative motion observed in numerical simulations .1^ 

As a future direction, it is promising to improve the 2D result by introducing 
nonlinear theory for Cd and Cr, for which Eq. (15.121) provides an encouraging prototype. 
On the side of SFD, more attention should be paid to cooperativity. In particular, 
it would be interesting to devise an extension of SFD in which a relaxation occurs 
cooperatively, as it would shed more light on studies of glassy dynamics than the 
classical SFD which is essentially /3 relaxation in ID systems. 
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